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We propose a method for quantization of Lagrangians for which the Hamiltonian, as a function of 
momentum, is a branched function with cusps. Appropriate boundary conditions, which we identify, 
insure unitary time evolution. In special cases a dual (canonical) transformation maps the problem 
into a problem of quantum mechanics on singular spatial manifolds, which we also develop. Several 
possible applications are indicated. 

PACS numbers: 



Quantum mechanics, like classical mechanics - and un- 
like, say, the standard model of fundamental interactions 
- is more a conceptual framework than a concrete account 
of specific physical phenomena. Physically interesting 
models based on quantum mechanics bring in additional 
structure. The analytical versions of classical mechan- 
ics have been a fruitful source of inspiration in this re- 
gard. There are procedures that allow one to pass from 
a wide variety of classical Hamiltonian systems to quan- 
tum models. In those constructions, guiding principles 
include Bohr's correspondence principle, which basically 
asks that the quantum model should reproduce, approx- 
imately, the original classical dynamics in appropriate 
limits; maintenance of symmetry; and internal consis- 
tency - specifically, unitary time evolution. There are 
also important but less fully developed, and possibly less 
rigorous, procedures using path integrals that allow one 
to pass directly from classical Lagrangian systems, in- 
cluding some singular ones, to quantum models In 
any case, we should regard the construction of quantum 
models as a creative process, open to innovation. 

Recently, in constructing models of possible time crys- 
tals Q , 0] , we were led to consider Lagrangians involving 
higher than quadratic powers of the time derivatives, and 
specifically the deceptively simple 

L = \i' \ i 2 (1) 

In the interesting case k > 0, the Hamiltonian for this 
Lagrangian is singular. Since the momentum involves a 
cubic in velocity, 

p = x 3 — kx (2) 

we can have either one or three real values of x corre- 
sponding to a given value of p, depending on whether 
\x\ — is positive or negative. Thus the energy func- 
tion 

*-£*-*-f**-jj*'. <3 > 



expressed in terms of p, is a multivalued function with 
cusps. See Figure 1. 




FIG. 1: Energy as a function of momentum. 

Here we propose methods for constructing quantum 
models corresponding to broad classes of classical sys- 
tems with branched structures in either momentum or 
position space. 

Branched quantization. Because the Hamiltonian 
based on Eqns. ([2] [3]) is not a single-valued function of 
p, and yet energy must surely qualify as an observable, 
p does not supply a complete set of commuting observ- 
ables. Therefore it will not be sufficient to label states 
with wave functions in (conventional) momentum space. 
On the other hand E is a single- valued function of x, so 
we might expect to construct wave functions ip(x). As x 
runs monotonically from — oo — > oo, p{x) evolves along 
the trajectory indicated in Figure 1, reversing direction at 
the cusps. This suggests that we consider wave functions 
that depend on p locally, but accommodate the back- 
tracking. Thus, denoting by p± = iy^f the cusp points, 
we have three components to the wave function, namely 
ipiip) for — oo < p ^ p+, ip2(p) for p~ ^ p ^ p+, and 
ip3 (p) for p_ ^ p < oo. Note that all three components 
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cover the range p- ^p ^ p+. 

A crucial issue is how the different branches join to- 
gether, i.e., what are appropriate boundary conditions. 
In this connection it is instructive to consider a more 
general class of Lagrangians than Eqn. {jTJ) , bringing in 
a quadratic potential V{x) — ^ax 2 . Directly from the 
Schrodinger equation we have an equation for the proba- 
bility density p^(p,t) = tp^ l (p,t)*tp fJi (p,t), defined on the 
branch labeled by p (where p = 1,2,3 for p_ < p < p + , 
p = 1 for p < and p = 3 for p > p+): 



dt 



ice, / + 5 2 ^ , v 



dp 2 



dp 2 
(4) 



using, in the second step, V{x) — > V(i-^). This sub- 
stitution implements the basic Heisenberg commutation 
relation, and also reflects the role of p as the generator of 
spatial translations. From Eqn. (j4]) we infer an equation 
of the current-conservation type for p = J2 Pn 



^P + ^l 
dt dp 



= 0; i - £t« 



ia , , * dtp. dip* 



dp dp 



O (5) 



Eqn. Q will lead to conservation of the integrated prob- 
ability J p if we can drop contributions from j at the 
endpoints. (Note that j receives contributions from two 
branches at each endpoint p±.) We also require our 
boundary conditions to be linear, so that our Hilbert 
space will support superposition, and so that they lead 
to a physically sensible eigenvalue problem for H . The 
choices 



4>2(p+) ; 



dip 



dp 



dtp 2 
dp 



(P+) (6) 



and their analogues at p_ manifestly give the required 
cancellation in j. Now consider the eigenvalue problem 
for the time-independent Schrodinger equation. For a ^ 
0, we get a second-order differential equation for tp(p). 
Thus, on each branch, for each value of energy, there are 
two disposable constants, making six altogether. Eqn. ([5]) 
gives us four constraints among these constants, and nor- 
malizability (absence of growing modes) at p — > ±oo 
gives us two more (one for each of ipx,ips). Thus the 
number of constraints matches the number of constants, 
as in conventional quantum potential theory. 

If a = the derivative conditions in Eqn. are 
not required and should not be imposed. On the other 
hand if V is a higher-order polynomial, we must require 
augmented boundary conditions. We will discuss those 
presently, after introducing a different (dual) viewpoint. 

Dual Viewpoint. One can hardly fail to notice that 
the manipulations we performed in momentum space, in 
connection with probability conservation, resemble ma- 
nipulations usually performed in position space. Thus it 
is natural to consider what our models look like after the 
substitution p — > x, x — ¥ —p, which preserves the struc- 
ture of quantum mechanics. After this substitution, our 
multi-valued kinetic energy becomes something perhaps 



less unconventional, that is, a multivalued potential. In- 
deed, we may think of a wire with kinks, as in Figure 
2a. Intermediate values of x are triply represented, and 
physical conditions will be different at different points 
along the wire, even if they are represented by the same 
x, so a branched wave function is manifestly appropriate 
to describe this physical system. 



(a) 



(b) 



(c) 



FIG. 2: (a) A wire with kinks, (b) A wire network with two 
junctions, (c) Box graph for a network with four junctions 
and a loop. 

From this dual point of view our quadratic potential 
V(x) — > V(p) = \ap 2 becomes the conventional ki- 
netic energy of a particle with mass m = 1/a, and the 
branched kinetic term becomes a multivalued potential 
W(x) in position space. Thus we have wave functions 
tpi(x) defined for —00 < x ^ x+ subject to W\(x), tp2{x) 
defined for X- ^ x ^ x+ subject to W^ix), and ^3(2;) 
defined for x_ ^ x < 00 subject to Ws(x), and boundary 
conditions similar to Eqn. (J6)), after the obvious substitu- 
tions of x for p. 

Now let us consider a quartic potential V(x) — x 4 + 
ax 3 + fix 2 +72;. In dual variables this leads to a kinetic 
energy that is a quartic polynomial in p, H = p A — ap 3 + 
ftp 2 — IP + W(x). We find a probability current (in the 
dual x space) 

I, ^d 3 tp dip^ d 2 tp d 2 ^ dtp d 3 ^ . 
^ i dx 3 dx dx 2 dx 2 dx dx 3 

, lf d 2 tp dtp'' dtp d 2 tp^ 

— otitp — + tp) 

dx 2 dx dx dx 2 



(7) 



where tp is a column vector with the requisite number of 
components tp^ in each momentum range. We will insure 
conservation of probability with the boundary conditions 



d n tp l 
dx n 



1 ' dx n 



< n < 3 



(8) 



at the kinks, provided that a and 7 change sign between 
the branches. This augmentation of the boundary con- 
ditions also leads to a good eigenvalue problem, since we 
have both twice as many disposable constants and twice 
as many conditions as in the quadratic case. 
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Inspired by the wire analogy it is natural, and for later 
purposes instructive, to consider networks analogous to 
the geometries of electric circuit theory, where we put 
quantum dynamics on graphs [4]. Let us consider what 
is required to insure no flow of probability into a node 
where several lines indexed by n come together. If the 
momentum dependence on each line is simply p 2 , and we 
orient each line so all coordinates x^ flow into the node, 
then the "Kirchoff" boundary conditions 



4>1 = 1p2 = ■■■ 



(9) 
(10) 



insure that no probability accumulates at the node. 
These natural conditions give good eigenvalue problems 
for the "Compton" tree graph and the box graph dis- 
played in Figure 2b, c (and many others). In the Comp- 
ton graph, Eqn. © yields 2x3 = 6 conditions at the two 
nodes, which together with four conditions at infinity for 
the external legs gives 10 conditions, as is appropriate 
to five lines with two disposable constants each. In the 
box graph, Eqn. (JSJ yields 4x3 = 12 conditions at the 
four nodes, which together with four conditions at infin- 
ity gives 16 conditions, as is appropriate for eight lines 
with two disposable constants each. Though it is simple 
and natural, this is by no means the only possible set-up 
consistent with the general requirements of the frame- 
work of quantum mechanics Q. We can vary both the 
Hamiltonians and the boundary conditions. 

To illustrate the possibilities, consider that each line 
might have its own kinetic Hamiltonian a^p 2 + (3^p (and 
of course its own potential). We can insure conservation 
of probability with boundary conditions of the form 



at the node, together with 

\ - * dip^ 
> a u K„ - = 



(11) 



(12) 



continuous. Similarly, the boundary conditions Eqn. (|SJ) 
for quartics unfold into continuity for ip{x) and its first 
three derivatives. 

Potentials. Now we return to our original problem, 
the issue of quantizing the kinetic Lagrangian of Eqn. (^Q) 
allowing for a general potential V(x). Inspired by the 
preceding unfolding procedure, we formulate our wave 
function in term of a variable that locally reduces to ±p 
plus a c-number, but covers all three branches following 
the same flow directions as d>. Thus we introduce 



p-p++p- 

-P + P++P- 
p + p + -p- 



for £ ^ p_ 

for p_ < £ < p + 

for p + ^ £ 



(14) 



and the decomposition of wave functions 



V(0(i-#(£-p-)) 
MO + MO + MO 



m-p+)) 



(15) 



where H is the Heaviside function. In this formulation p 
is realized (piecewise) as a modified multiplication opera- 
tor, with slightly different modifications on each branch. 
We can use that fact to write the operator V as an explicit 
kernel in £ space. Thus we transform tpi(0 to x space, 
where V acts as multiplication, and transform back as 
follows: 



u[x) = 



(Vu)(x) 

(VMO 



e lpx M0)^ 



V(x)u{x) 



27T 



V(x) dx 



(16) 



so long as /3 M | k m | 2 =0. With quartics, the possibilities 
proliferate. 

Unfolding. Quadratic models with branching in the 
dual viewpoint can be unfolded, upon which the eigen- 
value problem assumes a conventional form. In equa- 
tions: If we write define a real variable x with 



X 
X 
X 



x — x+ + X— for x ^ x — 

x + x+ + X- for X- x ^ x -\ 

x + x+ — X- for x+ ^ x 



(13) 



then as x evolves monotonically from — oo — > oo it covers 
each branch of x uniquely (as x covered p in our ear- 
lier discussion). The boundary conditions Eqn. ©, with 
p — > x become, upon transcription into the unfolding 
variable x, the statement that ip(x) an d its derivative are 



Note that the result of V acting on ipi generally does not 
vanish for £ > p_. It is realized as an operator of the 
Wiener- Hopf type. 

Performing the same manipulations on in the other 
intervals, we arrive at K-i = K*, K3 = K\ and 



(VMO = j (k^ - 0)M0) + k 2 (M 0)M0) 

+K 3 (£ - OM0)) d0 (17) 



The peculiarity of Ki arises from the reversed flow of p, as 
a function of £, in the medial interval. In the symmetric 
case V(x) = V(—x) all the Ks are real and equal, and V 
becomes an ordinary convolution operator. 

If V(x) is not symmetric, however, we must reconsider 
our procedure, because the V defined in Eqn. (II T[) is not 
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Hermitean. Indeed, although each Kj satisfies the her- 
miticity condition Kj{£' , £) = Kj the full kernel 

+ K 2 (Z',t)(H(£-p-)-H($-p + )) 
+ K 3 (Z',t)H(£-p + ) (18) 

does not. Thus to reach a consistent quantization we 
must impose K2 = K\ also (and not K2 = K*). The 
sign changes for a and 7 required in Eqn. ([8]) foreshad- 
owed this conclusion. By adopting this modified quanti- 
zation condition, we lose the Heaviside functions and ar- 
rive at a (manifestly Hermitean) convolution. The mod- 
ified quantization condition entails that the basic com- 
mutation relation [£, x] = —i involves the unfolded £, not 
the mechanical p. 

For long-range potentials V(x) the formal definition 
of V by Fourier transformation leads to derivatives of S 
functions, which must be defined through integration by 
parts on the momentum-space wave functions. In this 
way we make contact with our earlier discussion of poly- 
nomial potentials, and see why smoothness conditions 
connecting the different zones in ip are required, that be- 
come more demanding as the order of V increases. 

The eigenvalue problem for H, written in the £ repre- 
sentation, generally leads to an integral equation, rather 
than the boundary value problem for a differential equa- 
tion, because H is not polynomial in momenta. This can 
be solved directly. Alternatively, energy eigenstates i\)e 
with energy E are characterized by the conditions 

(il) B \[H,0]\^ E ) = 
(tP E m E ) = E 

{M^e) = 1 (19) 

for a complete basis of operators 0, or by the variance 
conditions 

(iP E \H 2 \ip E ) (ip E \OOt\ip E ) + (ip E \HOO^ H\ip E ) = 



p is invertible: 



(ip E \H\i> E ) = E 
(tp E \ip E ) = 1 



(20) 



We can solve these approximately by iterative (Newton) 
root-finding and minimization methods, respectively, as 
explained in detail in [j| . The required evaluations of ex- 
pectation values, given our explicit expressions, provide 
a practical approach to the eigenvalue problem. 

Eigenvalue problems for several examples of branched 
quantization with different potentials V(x) are analyzed 
in @. 

Another unfolding method. Another natural choice of 
an unfolding coordinate for p is x. The phase space co- 
ordinates (x, x) are noncanonical, but one can formulate 
a symplectic structure and a Poisson bracket for them 
QQ, which reduces to the standard Poisson bracket in 
the momentum range for which the map between x and 



3i 2 



8F_dG 

dx dx 



d£dG 
dx dx 



(21) 



With this Poisson bracket and the Hamiltonian 



H(x,x) = ^x 4 -|i 2 + V(x), (22) 

one obtains Hamilton's equation F = {F,H}, which re- 
produces the equation of motion derived from Eqn. (flj . 
This Hamiltonian formulation demonstrates that time 
evolution is well-defined (except for when x = 
where the symplectic structure degenerates and nonde- 
terministic motion may occur [3) and offers an alterna- 
tive approach to quantization [6( . 

Comments: 

1. A very common and fruitful procedure in analyz- 
ing quantum many-body problems, is to model the 
effect of interactions on a given particle by an effec- 
tive one-body Hamiltonian (or Lagrangian) , solving 
the one-body problem, and constructing a many- 
body wave function as a suitable product, e. g. a 
Slater determinant. By widening the class of can- 
didate one-body Hamiltonians, we can hope to ex- 
tend this sort of analysis to wider classes of systems. 
"Swallowtail" structures similar to Figure 1 have 
appeared in the description of Bose-Einstein con- 
densates in lattice traps Q. Dynamical mean field 
theory 0] generates complicated time dependence 
in a time-translation invariant energy functional, as 
a consequence of interactions. Polynomial trunca- 
tion of such time dependence, using substitutions 
of the type (x(t) - x(t - 8)) n -> 5 n x n 0, with re- 
tention of spatial structure, brings us to the sort of 
models considered here. 

2. A particularly interesting case arises for periodic 
potentials V{x). In that case, famously, conven- 
tional kinetic terms lead to band structures: The 
energy becomes a multivalued function of the quasi- 
momentum. Our branched Hamiltonian already 
for V = has a sort of band structure, associ- 
ated with the branches, in a region of momentum 
P- < p < p + where the limiting are deter- 
mined by the form of the kinetic energy, not by 
any spatial periodicity. With a periodic potential 
added, both sources of banding are effective. Espe- 
cially interesting is the possibility of describing dy- 
namically induced insulating behavior (Mott phe- 
nomenon) at filling fractions determined dynami- 
cally by the value of k. 

3. Our earlier work on time crystals was somewhat 
schizophrenic. In the classical case we found 
systems with motion in their ground state using 
kinetic Lagrangians of the type considered above. 
In the quantum case [|| , not knowing how to treat 
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such Lagrangians, we used a different mechanism, 
based on a more conventional kinetic term, that de- 
pended on the discreteness of (generalized) quan- 
tum angular momentum. A possible experimen- 
tal realization has been proposed [10(. With the 
method here described we can quantize the classi- 
cal time crystals, and construct much more general 
candidate models of quantum time crystals. 

4. We initially passed to the dual models to guide 
our intuition, but they appear to have consider- 
able independent interest. The central observa- 
tion, that appropriate, fairly simple boundary con- 
ditions, both on wave functions and (odd) inter- 
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actions, appears to allow construction of unitary 
quantum mechanics on a wide variety of singular 
manifolds. One can even relax the boundary con- 
ditions, at the cost of allowing probability to flow 
into and out of the designated points - in other 
words, by allowing the points to have internal de- 
grees of freedom. This procedure appears natural, 
specifically, in the modeling of black holes, where 
(in the Euclidean formalism) the horizon appears 
as a sphere attached to a point. 
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